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We present an exact calculation of the coherent thermal conductance in a 1-D multilayer photonic 
crystals (PC) using the S-matrix method. In particular, we study the thermal conductance in a 
bilayer structure of slabs of Si/vacuum or A^Oa/vacuum by means of the exact expression for the 
radiative heat flux. We compare our results with results obtained in previous works. Our results 
show that the coupling of surface modes as well as material losses play a fundamental role in the 
definition of the thermal conductance of PCs 
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Recently there has been a growing interest in exploring nanoscale heat transfer theoretically and experimentally 
which is triggered by the fact that radiative heat flux at the nanoscale can be much larger than that between two 
black bodies pj and quasi- monochromatic 0, 0| which makes it very promising for near- field thermophotovoltaics [IJ . 
The tremendous increase in the amount of transfered energy for distances much smaller than the thermal wavelength 
(Ath = frc/k^T) which can be several orders of magnitude larger than the value predicted by Stefan-Boltzmann's law 
can for dielectrics be attributed to the contribution of a large number of coupled surface phonon polariton modes 0, 0] ■ 

The contribution of surface modes (SMs) is indeed very important for nanoscale heat fluxes and many researchers 
have tried to enhance the amount of transfered heat by using this effect. Volokitin and Persson YA have pointed out 
that thin metallic coatings on a substrate can increase the nanoscale heat flux, Biehs et al. 8, 9] and Francoeur et 
al. [l(| have shown that one can use the cou plin g of SMs in thin metallic or dielectric films to enhance the nanoscale 
heat flux, Ben-Abdallah and coworkers [Tl| - |l4| and Francoeur et al. [H[ have also considered this effect between 
two finite slabs or media with several layers, Fu and Zhang [l6| have studied how doping affects the surface mode 
contribution, van Zwol et al. (l7j have shown that large nanoscale heat fluxes in phase change materials are due to 
SMs, and Svetovoy et al. [HI and Ilic et al. [l9| have pointed out that thin sheets of graphene allow do control or 
modulate the surface mode contribution. Very recent works have also considered heat fluxes for artificial structures 
and/or meta-materials supporting SMs in the infrared regime [20l - |2~^ |. or the surface mode coupling in many particle 
systems [25| . 

In this letter, we will revisit the theory of thermal conductance by photons within a PC as depicted in Fig. [TJ We 
will provide an exact expression for the thermal conductance inside a 1-D PC for arbitrary dispersive and dissipative 
material slabs. In particular, this allows us to determine the transmission coefficients (TC) for the Bloch states inside 
the PC. In previous works Lau and co-workers (2(| have assumed that the TC equals its maximum value of one when 
losses can be neglected. Surprisingly, by comparing the exact results of our calculation with the results of Ref. (2(| 
we find that in the limit of vanishing losses the TC for the total internal reflection modes goes to zero and not to 
its maximum value. In fact, we find that the TC is very sensitive to the losses inside the PC slabs. In addition, our 
exact expression takes the contribution of SMs to the thermal conductance inside the PC into account as well. We 
will show that this surface mode contribution can be crucial for the thermal conductance inside a PC. 
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Figure 1. Sketch of the situation considered here. Two identical 1-D bilayer PCs are separated by a vacuum gap of distance d. 
Here we choose d — h to determine the thermal conductance inside an infinite 1-D PC. 

In order to derive the expression for the radiative heat flux inside a PC, we assume that we have first two semi- 
infinite PCs as depicted in Fig. [T] The bilayer structure has a period of a = l-y + 1% where l\ is the thickness of the 
material layer with complex permittivity e\ — e[ + ie" and I2 is the thickness of the vacuum layer with 62 = 1. The 
heat flux 

^ r = M2i,d)Ar lr (1) 

between two such semi-infinite structures having a temperature difference A7j r = T r — T\ across the vacuum gap of 
distance d can be derived from the expression given by Polder and van Hove The heat transfer coefficient (HTC) 
within an infinite PC can then be obtained by setting d = I2 



^ f°°dw 96 (T) f 



1= s,p 



^Ti(w, K ;d=h). (2) 



Here, the time derivative of the Bose-Einstein function is given by dO(T)/dT = (huj) 2 / '(k B T 2 )e hu - '> lk ^ T / '(e^ '> ' ksT - l) 2 
and evaluated at the temperature T\ of the last slab of the PC at the left hand side. The TCs %(ui, k; d) for s- and 
p-polarized waves (i — s, p) are given by [l[ 

Ti{u,K]d) - { . , 1 ' (i jr )e -2|fc, |i (3) 
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where D- r = (1 — R\R\e 2lkzad )~ 1 is a Fabry-Perot-like denominator with k 2 — to 2 /c 2 — n 2 and n 2 — k 2 + k 2 . R s 
and R p are the reflection coefficients for the two semi-infinite PCs and can be calculated with the standard S-matrix 
method for layered media [HI EH Hz( ■ Note, that the TC is for propagating modes with parallel wave vectors k < uj/c 
different from the expression for evanescent modes with parallel wave vectors k > uj/c. 

Now, we are in a position to compare results from the exact expression in Eq. ([2]) with the results in Ref. [26j . First 
we note, that in the approach in Ref. [26[ the authors assume that the TCs equal their maximum value of one for all 
propagating bloch modes inside the PC. That means, the integral over all parallel wave vectors k is replaced by 

i—s,p 

Here the prime notes that the integral is for each frequency u> carried out over the whole parallel wave vector range 
which allows for propagating solutions inside the PC, i.e., over the photonic Bloch bands. The photonic Bloch bands 
can be determined from the dispersion relation for the Bloch modes (see Ref. Hence, the results from Refs. [26[ 

for the coherent thermal conductance give the upper limit for the contribution of the propagating Bloch modes. As 
we will see in the following, the exact result can be very different from such a calculation due to losses, resonant SMs 
as well as evanescent Bloch modes. 
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Figure 2. (Color online) (a) The TC for s-polarized modes defined in Eq. (0 in to-tt plane for ei = 12 + 0.001 ■ i. In magenta: 
numerical results for the boundaries of the Bloch mode dispersion relation [27]]. (b) A(u>) from Eq. Q normalized to A vlic — 
(lo I c) 2 / (2ir) 2 . The solid line represents the result for = 12 and 7s — T v = 1 for Bloch modes. The dashed lines show the 
exact results for fixed e[ — 12 and different e". (c) The HTC for the same permittivities as in (b) normalized to the black body 
result. 



In Fig. [UJa) the transmission coefficient T s is plotted in the lu-k plane choosing t\ = 12 + i ■ 0.001 for a PC with 100 
slabs. It can be seen that although the imaginary part of the permittivity is very small, corresponding to a system 
with vanishing losses, the TC is less than one for most parts of the Bloch bands. We find similar results for the 
p-polarized modes. In Fig. HJb) we present the numerical results for A(uj) when integrating the TC over k using Eq. 

The plotted values are normalized to the maximum value possible for propagating modes A vac = ^(ui/c) 2 inside 
the vacuum gap. The solid black line represents the result from Ref. [26j and the colored dashed curves represent the 
exact results using the same e[ as in Ref. [26[ but for different e'{. The best agreement with the black curve is found 
for e'l = 0.1. When decreasing the losses by making e'{ smaller than e± =0.1 then A(uj) decreases as well for nearly 
all frequencies so the deviation from the black curve gets larger. This means that for vanishing losses the TC does 
not converge to its maximum value for all Bloch modes. On the other hand, when making e" larger than e'{ = 0.1 the 
skin depth S s = l/[w/clm( v / ei)] inside the material slabs decreases and attains for e" = 10 (dashed magenta curve) 
values on the order of the period a of the PC so that the field is damped at this scale. Hence, the heat flux is not 
coherent anymore and the Bloch band structure in T s / P disappears. In fact, then the heat flux is due to Fabry-Perot 
modes of the cavity formed by separation gap and explaine the smooth and weak oscillating behaviour of A(u>) when 
e'l = 10 [Fig. Efb)]. Finally, in Fig. [2jc) we show the HTC h\ r versus the temperature T for the same permittivities 
as in Fig. [2fb) normalized to the vacuum or black body value h vac = (7r 2 fcgT 3 )/ (15c 2 H 3 ). Especially, for small T the 
deviation with respect to the results of Ref. (2(| is relatively large, whereas for large temperatures we obtain values 



very close to the 'universal' value found in Ref. [26 1 
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Figure 3. (Color online) (a) The TC for p-polarized modes defined in Eq. <[3j) in uj-k, plane for AbOa/vacuum structure with 
h = li = 5 /im. In magenta: numerical results for the boundaries of the Bloch mode dispersion relation [271 ]. The white dashed 
lines mark the frequency regions Ai and A2 where coupled SPPs exist, (b) HTC h\ r normalized to h[ r where A(uf) = A'{uj) 
and plotted vs. the vacuum gap 1%. For the solid curve the AI2O3 slab l\ varies analog to I2 and the dashed curves show the 
results for a fixed h = 5 /im and h — 50 nm while varying fa- 



Here below we examine the behaviour of structures which are able to support SMs, surface phonon polaritons 
(SPPs). To do that we consider an A^Os/vacuum PC with 100 slabs at 7} = 300 K. For this material combination 
SPPs not only exist for the p-polarized modes but also play the important role for heat transfer at subwavelength 
distances. In Fig. [3^a) we have plotted the TC for p-polarized modes in uj-k plane. It is obvious that not only the 
Bloch modes and Bloch SMs contribute to the heat conductance but also coupled SPP modes which can be identified 
in the frequency bands Ai and A2 where ei < — 1. To compare our exact calculations with results from (2(| for the 
AI2O3/ vacuum PC we have plotted in Fig.[3fb) the HTC h\ r versus the vacuum gap I2 for different l\. The results are 
normalized to the HTC h' h from Ref. [2fJ for which A(uj) = A'{lo) [see Eq. Q]. It can be seen that the exact HTC 
can be nearly four orders of magnitude larger than the HTC calculated with the approximative method at I2 — 10 nm. 
This can be attributed to the SPP mode contribution which is proportional to 1/1% for small gap sizes. 
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